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ABSTRACT 


This  study  was  undertaken  by  the  Research  Directorate,  Weapons 
Laboratory  at  Rock  Island,  to  develop  a  digital  computer  program  by 
which  the  two-dimensional  temperature  profile  in  gun  tubes  can  be  com¬ 
puted  under  realistic  physical  conditions.  A  mathematical  model  was 
presented  in  which  variable  geometry,  temperature-dependent  thermal 
properties,  and  variable  conditions  at  the  boundaries  are  considered 
A  numerical  algorithm,  In  which  the  method  of  explicit  finite-differ¬ 
ences  is  used,  was  developed  for  the  mathematical  model  and  was  pro¬ 
grammed  for  the  digital  computer.  A  numerical  example  was  computed  to 
check  the  computer  program.  The  program  and  all  subroutines  functioned 
properly.  No  numerical  Instability  nor  convergence  problems  were  en¬ 
countered. 
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INTRODUCTION 


Small  caliber  automatic  weapons  are  subjected  to  extremely  high 
operating  pressures  and  temperatures.  Energy  from  the  hot  propellant 
gas  is  absorbed  by  the  gun  tube  at  a  much  faster  rate  than  it  is  dissi¬ 
pated  to  the  surroundings.  Temperature  rises  occur  quite  rapidly,  and 
result  in  erosion  or  loss  of  strength  of  the  gun  tube  material.  High 
pressures  may  cause  the  tube  to  become  ruptured  when  the  temperatures 
are  increased  sufficiently.  The  Important  point  is,  therefore,  that 
gun  tube  designers  be  able  to  predict  the  temperature  distribution  of  a 
particular  gun  tube  design.  The  purpose  of  the  present  work  is  to  de¬ 
velop  a  digital  computer  program  by  which  gun  tube  temperatures  can  be 
computed  for  physically  realistic  conditions.  These  physical  conditions 
involve  variable  axial  geometry,  temperature-dependent  thermal  proper¬ 
ties,  variable  firing  schedules,  and  variable  thermal  boundary  condi¬ 
tions  at  both  the  bore  and  exterior  surfaces.  Results  obtained  from 
the  computer  analyses  may  be  used  to  determine  areas  of  excessive  tem¬ 
perature  rise,  to  estimate  maximum  burst  time,  to  provide  information 
necessary  for  thermal  stress  analyses,  and  to  indicate  necessary 
changes  for  improved  thermal  performance. 
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MATHEMATICAL  MODEL 


In  this  section,  the  physical  mechanisms  of  heat  transfer  from  the 
propellant  gas  to  the  gun  tube  are  described  and  a  mathematical  formul ac¬ 
tion  is  given  for  the  temperature  distribution  in  the  tube.  An  illus¬ 
tration  of  the  gun  tube  is  shown  in  Figure  1. 

The  gun  tube  material  is  considered  to  be  isotropic,  but  the' ther¬ 
mal  properties,  p(T),  C(T),  and  K(T),  are  known  functions, of  temperature. 
The  assumption  is  that  angular  temperature  variations  are  small,  com* 
pared  with  radial  and  axial  temperature  variations.  Thus,  only  a  two- 
dimensional  temperature  field  must  be  considered.  Heat  flows  from  the 
hot  propellant  gas,  whose  temperature  is  represented  by  Tg(r,z,t),  to 
the  tube,  whose  temperature  is  denoted  by  T(r,z,t).  The  assumption  in 
this  analysis  is  that  the  heat  flux  from  the  gas  to  the  bore  surface  is 
specified  or  that  the  propellant  gas  temperature  is  a  known  function  of 
time  and  position,  and  that  a  heat  transfer  coefficient,  hjCRj.z.tJ) 
exists  which  is  also  a  known  function.  And  finally,  the  assumption  is 
that  continuous  variations  in  the  outside  diameter  of  the  tube  may  be 
adequately  approximated  by  a  finite  number  of  step  changes  in  the  exte¬ 
rior  diameter.  The  actual  diameter  as  being  approximated  by  three  step 
changes  is  shown  in  Figure  1.  The  number  of  step  changes  may  be  greater 
or  less  than  three,  dependent  upon  the  situation.  Close  approximation 
of  any  taper  of  the  outside  diameter  by  use  of  a  greater  number  of  step 
changes  is  possible.  The  analysis  will  be  illustrated  with  the  use  of 
three  step  changes.  However,  the  computer  program  was  written  so  that 
any  desired  number  of  step  changes  in  the  external  diameter  could  be 
handled.  , 

The  governing  partial  differential  equation  for  the  gun  tube  is 
given  by 


3T  _  K(T  r 
it  '  p(T)C(T)  [ 


92T 


+  1 


£T 

3r 


92T, 


x  1  9K  r  /  3T\  •>  ,  /  3T\  , 

+  p(T)C(T)  3T  C(3F  +  (3T}  ]  ( 

Equation  1  is  nonlinear  due  to  the  presence  of  temperature-dependent 
properties.  The  boundary  conditions  for  the  tube  illustrated  in  Figure 
l  are 


r  =  Rj ,  0  :  Z  :  Z3  , 

-K(T)  |I  (Rj.Z.t)  =  q"(Ri,Z,t)  12J 
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ACTUAL  TUBE  PROFILE 


(3) 


r  =  R4,  0  4  Z  *  Zj 

-K(T)  ff  (R4,Z,t)  =  h(T)[T(R4,Z,t)  -  To] 

+  e(T)o[T‘*(R4,Z,t)  -  To**] 

r  ■  R3,  Zj  i  Z  *  Z2 

-  K(D  ff  <R3.Z.t)  =  h(T)[T(R3,Z,t)  -  To] 

+e(T)o[T4(R3,Z,t)  -  To4]  (4) 

r  =  R2,  Z2  -  Z  4  Z3 

-  K(T)  f£(R2,Z,t)  =  h(T)fT(R2,Z,t)  -  To] 

+e(T)o[T4(R2,Z,t)  -  To4]  (5) 

0  <  R  <  R4,  Z  =  0 

-  K(T)  |f  (r.O.t)  =  q"(r,0,t)  l6j 

R3  -  r  s  R4,  Z  =  Zj 

-  K(T)  |f  (r.Zj.t)  =  0  (7, 

R2  i  r  1  R3.  Z  =  Z2 

K(T)  r *Z2 ,t )  =  0  (gj 

Ri  -  r  £  R2,  Z  =  Z3 

-  K(T)  ff  (r,Z3,t)  =  h(T)[T(r,Z3,t)  -  To] 

+E(T)o[T4(r,Z3,t)  -  To4]  (9j 
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The  Initial  conditions  are 
T(r,Z,0)  -  T1  (r,Z) 


(101 


The  heat  flux  given  In  Equation  2  may  be  a  specified  function,  or 
may  be  expressed  In  terms  of  a  heat  transfer  coefficient  and  the  'ocai 
difference  between  the  bore  surface  temperature  and  the  bore  center 
line  gas  temperature.  In  the  latter  case,  the  boundary  condition  is 
given  by 


q"(Rj  ,Z,t)  ■  h(R1,Z,t)[Tg(O.Z,t)  -  T(R, ,Z.t)3 


(11) 


The  heat  flux  given  In  Equation  6  must  also  be  specified  Th<&  surface 
may  lose  energy  to  the  surroundings,  or  exchange  heat  with  some  oner 
section  of  the  weapon.  Boundary  conditions  for  this  surface  must  be 
specified  on  an  Individual  basis. 

The  temperature-dependent  thermal  properties  given  in  all  the  gov¬ 
erning  equations  must  be  evaluated  at  the  temperature  of  the  po'nt  at 
which  the  equations  are  being  evaluated.  The  heat  fluxes  at  the  bound¬ 
ary  locations  at  which  step  changes  have  been  used  to  approximate  con¬ 
tinuous  variations  In  external  diameter  are  assumed  to  oe  in  the  rad’ai 
direction  only.  This  Is  indicated  In  Equations  7  and  8.  The  relation 
form  factor  for  all  other  external  surfaces  has  been  taken  as  unty 

The  set  of  equations  given  above  cannot  be  solved  ana’yt  ta'V,  so 
numerical  techniques  must  be  employed.  The  method  of  explicit  finite 
differences  was  chosen  to  solve  the  equations.  The  details  ot  the  nu¬ 
merical  algorithm  are  given  In  the  following  section 
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NUMERICAL  ALGORITHM 


To  determine  the  temperature  distribution  In  the  gun  tube,  the 
tube  Is  first  subdivided  Into  a  finite  number  of  discrete  lumps.  The 
subdivision  of  the  gun  tube  Is  Illustrated  In  Figure  2.  The  tube  has 
been  divided  Into  three  sections  In  both  the  axial  and  the  radial  di¬ 
rections.  The  number  of  nodes  In  the  first  radial  section  Is  ij.  In¬ 
cluding  the  Interface  between  sections  1  and  2.  The  second  and  third 
radial  sections  contain  12  and  13  nodes,  respectively.  In  the  first 
axial  section,  j,  nodes  are  present.  Including  the  Interface  node. 

The  second  and  the  third  axial  sections  contain  j2  and  j3  nodes,  re¬ 
spectively.  The  spatial  Increments  between  the  nodes  are  given  by 


Af,  ■ 

(R2-Rj) 

1TFTT 

(12) 

ar2  ■ 

(R3-R2) 

~TT~ 

(13) 

ar3  « 

(Rv-R3) 

(141 

aZi  ■ 

h 

(15) 

TYFR 

aZ2  ■ 

(Z2-Z.) 

J2 

(16) 

aZj  ■ 

Urh) 

Jj 

(17) 

The  nodal  point  locations  are 

rn  ■  K>  1  s"5  1. 

(n-lj) 

rn  *  Rz  ♦  X2  (Rj-M.  1i  <  n  s  1j+12 
(n-1 )-12) 

rn  .  r3  ♦  - j- - (R„-R3),  V12  *  n  i  1i+12«13 


tlfil 

(19) 

(20) 
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FINITE  DIFFERENCE  GRID  NETWORK 


1 


1 


(m-ji) 

zm  =  zi  +  -3^—  (Z2-Z1).  Ji  <  m  -  J'i+32 

zm  =  Z2  +  - 33 -  (Z3'Z2).  Jl+J2  <  m  -  Jl+J2+J3 


(21) 

(22) 

(23) 


1  ' 

The  temperature  of  each  lump  Is  assumed  to  be  uniform  and  equal  to 
the  temperature  of  its  center  or  nodal  point.  The  spatial  derivatives 
appearing  in  the  governing  equations  are  approximated  by  finite  differ¬ 
ence  relations,  which  have  been  determined  from  the  simultaneous  solu¬ 
tions  of  truncated  Taylor  series  expansions.  Equation  1,  written  in 
discrete  difference  notation,  is  -  1 


(ft) 


n,m 


3t'  o  C 
n,m  ^n,m  n,m 


m  +  r  (£>  ] 

n,m  n  n,m 


*  Z — T —  [(i7>2  *  (a7,!  ]  ■ 

pn,m  Ln,m  dr  n,m  V  n,m 


(24) 


Subscript  n  denotes  the  radial  node  location  and  subscript  m  denotes 
the  axial  node  locations.  The  subscripts  of  the  property  values  denqte 
that  they  are  evaluated  at  the  temperature  of  node  r,m. 


K 


n,m 


K  (T 


n,m 


(26) 


1 

The  boundary  conditions  must  also  be  written  in  difference  form, 
example.  Equation  2  is  written  as  1  1 


-K  (—1 
1 ,m  'ar'i.m 


1  ,m 


For 


(26) 


The  spatial  derivatives  In  Equations  24  and  26  are  approximated  by 
finite  difference  relations.  A  set  of  explicit  algebraic  equations1  re¬ 
sult  for  the  time  rate  of  temperature  change  at  each  node.  These  rates 
of  change  are  multiplied  by  a  finite  time  Increment  to  find  the  temper¬ 
ature  at  one  increment  of  time  later.  The  procedure  is  repeated  until 
the  final  time  period  of  Interest  has  been  reached.  The  finite  differ¬ 
ence  expressions  used  to  approximate  the  spatial  derivatives  are  shown 
in  Tables  I  and  II.  The  particular  derivatives  shown  are  for  the 
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1 


I 


I 


.1 


radial  direction.  The  expressions)  for  the  axial  direction  are  of  the 
same  form,  where  the  spatial  increment  is  aZ,  the  subscript  n  is'  fixed, 
and  the  subscript  m  is  a  variable.  The  appropriate  expression  to1 use 
depends  on  the  node  of  interest.  For  interior  nodes.  Number  2  'n  Table 
I  and  Number  3  in  Table  II  are  used.  At  the  bore  surface.  Number  i  dr 
Table  II  is  used  for  the  second  derivative;  the  first  derivative  is  de¬ 
termined  from  the  boundary  condition,  The  correct  expressions  for  tee 
external  boundary  nodes  are  those  of  theteeat  flux  boundary  condition:, 
and  Number  4  of  Table  II.  The  derivatives  for  the  nodes  adjacent  to  the 
boundaries  are  given  by  Numbers  2  and  3  of  Table  1  and  Number  2  ot  Table 
II.  ,At  the  interface  nodes,  by  which  sections  are  separated  in  white 
the  increments  between  nodes  may  become  changed  in  size,  the  corre.t  ex-1 
,  pressions  are  Numbers  4  and  5  of  Table  I  and  Table  ll ,  respectively 

v  (  I 

1  The  maximum  time  increment,  by  which  numerical  stability  is  ensued 
for  the  linear  diffusion  equation  with  this  algorithm,  is  g’ven  by 


i 

where  aX  is  the  distance  between  nodal  i points  and  w  is 


At  -  cC 


Mi 

Kw  i 


w  -  4 


hAX  1 
K 


I 


I 


The  present  set  of  equations  are  nonlinear  and  the  maximum  aU'owab'e  t  •  me 
interval  would  be  expected  to  be  less  than  that  gilV»n  by  Equafon  27 

•  1  i 

The  sufficient  condition  for  numerical  stability  and  con-erjen-e  is 
that  both  the  first  and  the  second  laws  of  thermodynamics  be  f.at"  "ea 
The  satisfaction  of  the  first  law  was  verified  by  the  perf unmani,?  o-'  of,' 
energy  balance  on  the  tub6  after  each  time  interval.  No  attempt  Was  made 
to  check  the  satisfaction  of  the  second  law  The  belief  was  that  the 
satisfaction  of  the  first  law  provided, an  adequate  check  ot  numer’.at 
stabi lity  and  convergence.  . 

I  1  '  :  ’ 

The  explicit  ti.nite  difference  algorithm,  described  in  th  s  sect  on', 
was  programmed  for  the  digital  computer  The  computer  prog- am  s  ot- 
sc'ibed  in  the1 f o 1  lowing  sectiqn,  and  a  program  listing  is  y’-en  m 
Appendix  A. 
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FINITE  DIFFERENCE  EXPRESSIONS  FOR  FIRST  DERIVATIVES 
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TABLE  II 


DESCRIPTION  OF  COMPUTER  PROGRAM 


A  digital  computer  program  was  written  for  the  evaluation  of  the 
numerical  algorithm  described  in  the  previous  section.  The  computer 
program  comprises  a  main  program  and  eight  subroutines.  The  main  pro¬ 
gram  contains  the  input,  the  output,  the  logic  operations,  and  the 
computation  operations  for  temperature  changes.  Detailed  calculations 
are  performed  in  the  subroutines.  The  name  and  the  purpose  of  each 
subroutine  is  given  below: 

1.  CONV  -  Provides  external  convection  coefficients  and 

emissivities . 

2.  QZSUB  -  Contains  operations  to  compute  the  axial  heat 
fluxes  at  the  external  surfaces  due  to  both 
radiation  and  convection. 

3.  QRSUB  -  Contains  the  operations  necessary  to  compute  the 
radial  heat  fluxes  at  the  surfaces  due  to  both 
radiation  and  convection. 

4.  AXIDER  -  Specifies  the  operations  for  the  computation  of 
the  spatial  derivatives  in  the  axial  direction. 

5.  RADDER  -  Provides  computations  for  the  spatial  derivatives 
in  the  radial  direction. 

6.  DKDT  -  Gives  derivative  of  thermal  conductivity  with 

respect  to  time. 

7.  XKKS  -  Gives  the  calculation  of  thermal  conductivity  as 

a  function  of  temperature. 

8.  LINEAR  -  Gives  specific  heat  as  a  function  of  temperature. 

The  input  to  the  digital  computer  program  consists  of  seven  READ 

ments  whose  required  input  data  are  listed  below: 

1.  M  -  number  of  radial  segments 

N  -  number  of  axial  segments 

SIGMA  -  radiation  coefficient 

TS  -  ambient  temperature 

TTIME  -  termination  time 

-  iteration  number  at  which  printout  is  desired 

I)  -  number  of  radial  nodes  in  each  sepent 

-  number  of  axial  nodes  in  each  sepent 
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4.  RS(J) 

5.  ZS(I) 

6.  KRA(I) 

7.  X(I),  Y(I)  - 


radial  boundaries  of  each  segment 

axial  boundaries  of  each  segment 

number  of  radial  segments  in  each  axial  segment 

temperature  versus  specific  heat  data 


Thermal  property  data  in  the  subroutines  are  for  SAE  4130  steel. 

If  a  different  barrel  material  is  to  be  analyzed,  the  functional  rela¬ 
tionships  in  these  subroutines  must  be  changed.  The  emissivities  and 
external  convection  coefficients  in  the  present  subroutines  are  con¬ 
stant.  If  these  values  are  not  constant  for  any  case  being  investi¬ 
gated,  the  proper  functional  relations  must  be  added  to  the  subroutines. 


The  output  from  the  program  consists  of  the  time  and  the  tempera¬ 
ture  at  each  node  for  those  iterations  for  which  printout  is  desired. 

A  complete  listing  of  the  digital  computer  program,  along  with  typical 
input  and  output  data,  is  given  in  Appendix  A. 
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NUMERICAL  EXAMPLE 

A  numerical  example  was  computed  to  check  the  digital  computer 
program.  The  sole  purpose  of  computing  the  numerical  example  was  to 
ensure  that  the  program  and  the  subroutines  were  functioning  properly. 
No  specific  weapon  was  considered.  With  reference  to  Figure  1?  the 
geometric  dimensions  used  in  the  example  were 


Rj  =  0.625  inch 

(29) 

R2  =  0.845  inch 

(30) 

R3  =  0.940  inch 

(31) 

Ru  =  1.088  inch 

(32) 

Zj  =  3.16  inch 

(33) 

Z2  =  12.5  inch 

(34) 

Z3  =  42.0  inch 

(35) 

Thermal  property  data  for  SAE  4130  steel  were 
2.013,  of  Aerospace  Structural  Metals  Handbook.2 
to  the  curves,  which  are  given  below. 

obtained  from  Figure 
The  data  were  adjusted 

K  =  (28.3  -  0.0087T)  T  *  1420°F 

(36) 

and 

K  =  (10.39  +  0.00347T)  p^r  sy,  T  >  1420°F 

(37) 

The  density  variations  for  SAE  4130  steel  are 
lowing  mean  value  of  density  was  used. 

small ,  and  the  fol- 

p  =  490  lb/ft3 

(38) 
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Tabular  data  for  specific  heat  were  used  in  conjunction  with  a 
linear  interpolation  subroutine.  The  specific  heat  data  from  Figure 
2.015  of  Aerospace  Structural  Metals  Handbook2  are  given  below: 

T,  °F  C,  BTlI/lb  °F 


0 

0.108 

200 

0.112 

400 

0.125 

600 

0.132 

800 

0.150 

1000 

0.160 

1200 

0.185 

1600 

0.180 

2000 

0.180 

2200 

0.150 

The  temperature  of  the  surroundings  was  constant  and  equal  to 
70°F.  A  constant  convection  coefficient  of  5  BTU/hr  ft  CF  and  a  con¬ 
stant  emissivity  of  0.5  were  prescribed  at  the  external  boundaries.  An 
effective  mean  propellant  gas  temperature  of  2000°F  and  an  effective 
mean  heat  transfer  coefficient  at  the  bore  surface  of  200  BTU/hr  ft  F 
were  used  in  the  calculations. 

The  computer  program  was  run,  with  the  use  of  the  data  given  above, 
for  a  continuous  firing  burst  of  12.5  seconds.  All  portions  of  the  mam 
program  and  its  subroutines  functioned  properly.  No  numerical  insta¬ 
bility  nor  convergence  problems  were  encountered.  The  temperatures  at 
all  nodes  were  printed  at  time  intervals  of  approximately  one  second. 

The  computed  bore  surface  temperatures,  as  functions  of  time  and  posi¬ 
tion,  are  shown  graphically  in  Figure  3  The  bore  surface  temperature 
rises  more  rapidly  at  locations  in  which  the  barrel  wall  is  thinnest, 
as  would  be  expected.  The  effects  of  axial  temperature  gradients  are 
minor,  except  where  an  abrupt  change  occurs  in  the  external  diameter. 
This  indicates  that,  for  this  specific  example,  a  less  complicated  and 
less  expensive  one-dimensional  numerical  program  could  be  used  over 
most  of  the  axial  length  without  the  introduction  of  any  major  errors. 
The  two-dimensional  program  could  still  be  used  in  regions  where  a  step 
change  exists  in  diameter. 

A  summary  of  this  investigation  is  given  in  the  following  section. 
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FIGURE  3  BORE  TEMFERATURE  VERSUS  DISTANCE  FROM  BREECH 

'  -AS  A  FUNCTION  OF-TIME 


T 


I 


SUMMARY 

>  •  t  ! 

A  mathematical  model  to  determine  the  two-dimensional  temperature 
profile  in  a  gun  tube,  under  realistic  physical  conditions,  is  pre¬ 
sented.  Variable  geometry ,!  temperature  dependent  thermal  properties, 
and  variable  conditions  at  the  boundaries  were  considered  in  the  mathe¬ 
matical  model.  A  numerical  algorithm  was  developed  for  the  mathematical 
model  by  use  of  the  method  of  explicit  finite  differences.  The  numerical 
algorithm  was  programmed  for  evaluation  by  the  digital  computer.  A 
i  numerical  example  was  computed  to  check  the  computer  program.  The  pro¬ 
gram  and  all  its  subroutines  fgnctioned  properly.  No  numerical  insta¬ 
bility  nor  convergence  problems  were  encountered. 
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C  *  OICIIAk  CDMUIII  HOWM  TO  OtTHMIM  T M>-0 1 MN1I0NAL  T(MM- 

C  tlAIUM  MOMli  IN  CUM  TUMI  •»  TMf  MINOO  00  (IMLICIT 

C  FINIIt  DIFFfMKil. 

C 

DINfNlIQN  11*0,  Mi.  0*1*0.10),  011*0.101,  0*111.  DII1I,  MCI*0l. 

1  tCUOl,  01*01,  CMt*OI.  0141  *01  •  Ml  11*0).  1X1*01.  IM*t«OI, 

2  or iioi.  01*1*0.10),  01*21*0,101.  otmo.ioi,  01121*0,101, 

>  jjiio).  mm.  oiioi.  mu,  us<*i,  usioi,  jinn,  ***i*i, 

«  *1201,  11201 

COMMON  /ILK1/JSUNI0I,  M,  NOONT ,  NO,  I1UM.  Ill,  M,  SICMA,  TS 
*1*0  1.  M,  N.  1ICM*.  It.  ItINt.  MM.  m 
M*IMI  201,  M,  M.  II  CM*.  11.  IIINf,  MM 

201  FO*M»H*M  N  II,  1*.  «M  M  •,  II,  II,  ON  II CM*  ••  (11.0.  II. 

I  IN  IS  .,  1)2.1,  1*.  OM  IIIMC  -,  (12.1,  tl.M  MM  ■,  101 

C  MM  -  MMICN  |I(**I|CM  Mtl*(  001**11  V(  M«IMf  1. 

C  N*  •  M00UCU1  M«  IMOICAf (1  MMICN  IHBAIIONf  TIMM  MINI  0UI1  Ml  DftlMO. 

*(*0  21, I JJtl 1 1,  I  •  I,  Ml 
M*IMI  202.  IJJUII.  I  ■  I. Ml 

202  FOONAt I «M  JJl/l II  111  I 
•(*0  21.IU1IJI.  J  ■  I.  Ml 
MAIM!  201.  IllllJI.  J  -I. Ml 

201  FO*M*l I4M  LLt/l II  111  I 
MM|  .  M  «  | 

*1*0  2,  I  *1121.  J  •  I.  MM 1 1 
M*INI  20*.  lailJI.  J  ■  I.MMII 
20*  FOOMAI I  IN  •}/  I  It  11. Til 

C  AIL  tf QMfMtl  ■  1CCMCMI1  00 1  AA0UILT  HUtl  M**f  1*W 

•(*0  2.  I  fllll.  I  -  I.  Ml 
M*IMf  201,  llllll,  I  •  I, Ml 
201  FJ*N*II«M  11/  IKIl.ril 

C  All  ItCMtNH  *  ItCXMll  OUT  *11  ALL*  *0*1  NAM  SAX 

*(A0  21.1***111.  I  •  I.  Ml 
MAIN!  20*.  IK*AIII,  I  «  I, Ml 
20*  FOOMAI  |*M  *«*/  lllllll 

MAO  2.  I  Kill.  VIII.  I  •  1,101 

1  FOAMATI2IIO,  fll.l,  2*10.0,21 101 

2  FOOMAI I IF  10.0 1 
21  FOOMAI 101 101 

C 

C  M*|MI  HO.MAIII.  AAAI/I.  AAAI  II 

C  110  FO*MAII«M  «*AIII  ••IA.I01.0M  1*0121  ••  lA.IOl.  OM  ***111  ••  I A I 

C 

KIMI  •  M  ,  M  .  I 

C  *  Hill  II  1(11  IMAM  0*  (OVAL  10  N 
C  Null  MAVt  M  MUMMAOM  ISUqlll't. 

C  MAI  1 1  It  MUMM*  0*  0*0111  UCIIOMt  M>  (KN  0*1*1  MCI  I  ON. 

00  1  *• I ,M 


*001*1  0IMN1I0M. 

0*1*1  DIM  Ml  ION. 
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RRAD  >  RRAIRI 
JSUMtRI  «  0 
00  i  J  •  I.RRAD 
J  JSUNIRI  ■  JJSIJI  ♦  JSUNIRI 

c 

MINI  ))2. JSUMI II,  JSUNI2I,  JSUMI  Jl 

]»  FORMAT  1 10H  JSUMI l I  -.I4.10I.10H  JSUMI})  ■.l4.t01.IOM  JSUNI1)  -,141 

c 

I SUM  •  0 
00  4  I  •  I.  N 

4  t Sum  •  i sum  »  ust  n 

C  I  SUM  REPRESENTS  THE  TOMl  NUMBER  Of  *1141  NOOES 

C  JJS-  NUMBER  OF  R40UI  NOOES  IN  THE  RESPECTIVE  SEGMENTS. 

C  US-  NUMBER  OF  *11*1  NOOES  IN  THE  RESPECTIVE  SEGMENTS. 

C  N  -  NUMBER  OF  *11*1  SEGMENTS. 

C  N  -  MAI | MUM  NUMBER  OF  R*OIAl  SEGMENTS. 

C  Ttl.JI  -  TEMPERATURES. 

C  ORIIt  -  RtOIAl  INCREMENT  CHANGE  I  OR 1 1 1  <  DRIJI  <  ORIS)  I 

C  Of  1 1 1  -  *11*1  INCREMENT  CHANGE  ISTART  Ml  TH  D(l|l  •  Oil})  ■  OIIS)  I 

C  TRIM.  AND  ERROR  FOR  CORRECT  VALUES  FOR  OIIII'S. 

C  MS  I J I  -  BOUNDARY  RAOI I  FROM  BORE  TO  OUTSIDE. 

C  ISIII  -  -  A1IM  BOUNDARIES. 

C  IS  -  AMBIENT  TEMPERATURE. 

C  HGIITI  -  WHERE  IT  IS  TEMP  SUBSCRIPT  A  NO  HC  IS  THE  CONVECTION  COEFFICIEN 

C  T  OF  THE  GAS  COMPUTED  AS  A  FUNCTION  OF  TEMP. 

C  II  -  THERMAL  CONDUCTIVITY  AS  A  FUNCTION  OF  TENP. 

C  CPI  II I  -  SPECIFIC  HEAT  AS  A  FUNCTION  OF  TENP. 

C  RHOIITI  -  OENSITY  AS  A  FUNCTION  OF  TENP. 

C  CNISIIT)  -  EMISSIVITV  AS  A  FUNCTION  OF  TENP. 

C  I  MCI  I T  I  -  CONVECTION  COEFFICIENT  AS  A  FUNCTION  OF  TEMP. 

C  1HRIITI  -  MAOIATION  AS  A  FUNCTION  OF  TEMP. 

C  TIIME  IS  TINE  OF  TERMINATION 

C 

C  THE  FOLLOWING  COMPUTES  THE  AIIAl  ANO  RAOIAl  CHANGES  IOEITASI 

00  4  J  •  I,  N 
A  IJSIJI  ■  JJSIJI 
00  T  I  •  I.  N 
T  USUI  •  USIII 

Dllll  ■  {Sill  /  I  XL  SI 1 1  •  1.01 

00  I  l>  2.  N 

B  Dllll  •  I  (SI  1 1  -  (SI  1-111/  USIII 
C  NUT  COMPUTE  THE  RAOI  I 

(Will  •  IRSI2 1  -  RSIIII  /  I X  J  Sill  -  1.0) 

Rill  ■  RSI  1 1 
J1K  ■  2 
JJJ  •  0 

00  II  Jl-  I.  M 
IF  I  JR  .10.  II  SO  10  T 
OR  I  Jl  I  •  IRSIJ1MI  -  <1  SI  JX 1 1  /  IJSIJI  I 
JU  «  JRR  «  JJS  I  JR- 1 1 
«  JJJ  •  JJJ  »  JJS I  JR  I 
00  10  J  •  JRR.  JJJ 

10  RIJI  «  R I  J—  1 1  *  OR  I  JR  t 
J-  J  *  I 

RIJI  «  RSI  JR  •  II 

11  CONTINUE 

RIJI  ■  RSINMI 
TINE  •  0.0 
NRONT  ■  I 
ISAB  •  TS  ♦  4*0.0 

NN  •  0  • 
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.mm  -  jjsui  »  i 

MM l  •  N-l 
00  14  J  ■  2.NN1 
14  JIIIJI  •  JIllJ-ll  ♦  JJSUI 
IK  ■  N  -  LISINI 
IK  *  1 
NN  •  USUI 
00  I*  II  •  I.N 
JSU  •  JSUMI III 
'  00  17  I  ■  IK.NN 
00  17  J  •  1,  JSU 

17  ru.ji  ■  is 

IHII  .£0.  Ml  CO  TO  1C 
IK  •  IK  ♦  USIII) 

NN  a  NN  t  IISIII»II 
1C  CONTINUE 
M  CONTINUE 

c  1 

C  NEKT  CONFUTE  THE  MAXIMUM  TINE  INTERVAL. 

C  OK  ■  Mill  I  SINCE  Mill  <  Mill  M  Oltll  I 

OX  •  Ml!) 

C 

C  NON  LOCATE  MAXIMUM  XKIITI.  AMO  MINIMUM' XHI I Tl 

C  XKNAX  •  XKIITIMAX 

C  XHNINHa  XHI  IT  ININ 

KKMAX  ■  21.) 

XKNIN  ■  12. • 

XHNAX  a  240.0  1 

XHMIN  ■  207.0 

c 

XMA  a  4.  ♦  4.  •  XHNAX  7  XKNIN  •  OX 
C  NON  LOCATE  RHOIITI  MINIMUM.  CXIITI  MINIMUM. 

C  XHMIN  ■  XNOIITINIM 

C  CXMIN  -  CXIITININ 
XHMIN  a  4*0.0 
CXMIN  a  .10* 

c  , 

071 1 1  a  XHMIN  •  CXMIN  7  KMX  7  XKNAX  •  OX  •  OX 
XNR  a  4.  ♦  4.  •  XHMIN  7  XKNAX  •  I  MSI  2 1  -  MSI  1 1  1  7  IXLSIll  -  1.0) 

0TI2I  a  XHMIN  •  CXMIN  7  XNR  /  XKNAX  •  IRSIIl-XSI  II  l»»2  7  XLSI 1 1 

l  ••  2 

00  I*  J  a  2.  N 

XNR  a  4.  «  4.  •  XHMIN  7  XKNAX  •  I  MSI J* 1 1  -  MSI J 1 1  7  XLSI J I 
1*  DTIJ*|)  a  X  HON  IN  •  CXMIN  7  AMR  7  XKNAX  •  IRSIJMI  -  RSIJII  ••  2  / 

1  XLSIJI  ••  2 

XN2  a  4.  ♦  4.  •  XHMIN  7  XKNAX  •  tSIll  7  IXJSIH  -  1.0  I 

OHN»n  a  XHMIN  •  CXMIN  7  XN 2  7  UMAX  •  tSIll  ••  2  7  XJSI1I  •*  2 

00  1*0  I  a  2.  N 

XNt  a  4.  *  4.  •  XHMIN  /  XKNAX  •  ItSIII  •  IStl-lll  7  XJSII) 

1*0  OTII.NI  ■  XHMIN  •  CXM|N  7  XHI  7  XKNAX  •  USIII  -  2SII-III  *»  2 
1  7  XJSII)  **2 

C  NOW  MUST  DETERMINE  SMALLEST  OT  ANO  SAID  VALUE  HILL  OE  THE  TINE 

C  INCREMENT. 

IOT  a  o 

>0  IOT  a  IOT  »l 
I XT  a  IOT  *  1 

IXIOTIIOTI  .ST.  OTIIXTII  SO  TO  SO 
OTIIXT)  a  OTIiOTI 
SO  MINT  1  a  HINT  •  1 

IXI1XT  .Nf.  KINT1I  CO  TO  20 
C  HINT  WIST  OS  READ  IN  AS  IN  ♦  N  ♦  11 
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i  Of H  •  OTIKINT1I 

OTH  •  OTH  /  1.5  ' 

.TIME  «  TIME  ♦  OTH  •  3600.0  I 

C  NON  COMMUTE  ALL  VALUES  DEPENDENT  ON  TEMPERATURE 

C  SUBSCRIPT  IT  DESIGNATES  WHICH  TIME  INTERVAL 

C  XHCII.JI  MUST  EITHER  BE  A  TABLE  READ  IN  OR  A  FUNCTION  COMPUTATION. 

C  i  NEAT  COMPUTE  THE  HEAT  FLUXES. 

LL1  •  LLSI1I 

CALL  QRSUBI T •  OR. TSAB. JM.LL  SI 
C  •  ••••*'• 

CALL  01 SUB I TSAB. JJ$. LLS, T.QE .KRAI  1 

C  NOW  CHECK  RAOIAL  TEMP  AGAINST  AXIAL  TEMP  01FF  A  NO  SHOULD  BE  CLOSE. 

C  ONLY  NEEO  TO  CHECK  BORE  TEMPS  AXIALLY  AGAINST  ADJACENT  RAOIAL  TEMPS. 

CAU  RADDER(T,DR»0TR.0TR2*QR»JJStLLS) 

C  NEXT  COMPUTE  AXIAL  OEAIVATIVESJ:  ' 

CALL  AXIOEAIT.D2.DT2.OT22.Q2.LLSI 
C  COMPUTE  TEMP  CHANGE  AT  EACH  NQOE 

IFINXONT  .EQ.  II  GO  TO  AM 
IFIMODINKONT.NRI  .NE.  0)  GO,  TO  SIO 
C 

AM  PRINT  |  500,  TIME.  NkONT 

500  FUK1ATI//5X,  BH  TINE  •  .  E12.A,  10X,  |BH  ITERATION  NUMBER  .  15/1 
510  IFF  •  0 

DO  503  IK  ■ 1.  N  ' 

I  IN  «  IFF  ♦  f 

IFF  «  IFF  ♦  LLSIIK) 

JSU  •  J  SUM  II K  I 

00  ASA  J  -  I.  JSU  >  ' 

00  ASA  I  •  I  IN.  IFF 
TT  •  T  (  I  6  J  > 

CALL  OKOTITT.  XXK,  OKT.  XRHO.  XCP.X,  VI 

c  okdt  is  f hf  thermal  conductivity  ,  density,  c  specific  heat  subroutin 

'  |  I F (NKONT  .NE.  0)  GO  TO  BOS 

PRINT  002.  XXK,  XRHO,  XCP,  0TR2II.JI.  RIJI,  DTRII.JI,  DTE2II.JI, 

I  OKT,  DTE  1 1 , J I 

S02  FORMATIAH  XXK  «,  E15.T,  5X,  7H  XRHO  •,  E15.T.  5X.  AH  XCP  >.  E15.T, 

I  5X,  I2H  0TR2II.JI  -,  E 15. T  /  7H  RIJI  •,  EI5.7,  5X,  UH  DTRII.JI 
2*.  EI5.7.5X.I2H  0TE2II.JJ  E15.7,  5X.  AN  OKT  E15.7  /  5X,  12H 
1  S  DTE  1 1, J  I  «,  E 15. 7 1 

BOS  OTOT  -  XXK  /  XRHO  /  XCP  •  I0TR2II.JI  »  I.  /  IRIJII  •  OTRII, 

1  Jl  •  DTE2I  #.  J I  I  ♦  1.  /  XRHO  /  XCP  •  DKT  •  IOTRII.JI  ••  2  ♦ 

2  OTEII.JI  ••  21 

>  ASA  TlT.JI  •  TII.JI  »  DTOT  «  OTH 

IFINKONT  .EQ.  II  GO  TO  AAI  , 

IFIMODINKONT.NRI  .NE.  01  GO  TO  50S 
C 

I  AAI  CONTINUE 

DO  501  I  «  IIN.  IFF 
PRINT  502.  I 

502  FORMAT  1 5X.  1BH  AXIAL  LOCATION  •  ,  151 
PRINT  505.  ITII.JI.  J  ■  I.JSUI 

501  CONTINUE 
501  CONTINUE 

505  FORMAT  I  SI SX, E 12. Al I 
5A0  NKONT  ■  .NKONT  »  l 
C  IFINKONT  .LE.  5001  GO  TO  AO 
IFITIME  .LT.  TTIMCI  GOTO  AO 
7 AO  CALL  EXIT 
ENG 
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SUBROUTINE  CONVI TT, XXHC, EMISS, XK ) 

EMISS  ■  .5 

XXHC  •  5.0 

CALL  XKKSITT.XK.DKTI 

RETURN 

END 


SUBROUT  INE  Q2SUBI TSA6, JJS.LLS, T,Q2,KRA I 

COMMON  /BLXI/jSUMIBIf  N,  NKONT,  NR.  ISUM,  LL1,  M,  SICMA,  TS 
DIMENSION  02140.30).  TI40.30I,  JJSI8I,  LLSIS).  KRAIS).  XI20). 

1  VI 20 1 
00  90  KR«1,N 
ISU  *  0 

DO  5  IK  •  I.  KR 
5  ISU  •  ISU  *  IASI  IX) 

IFIKR  .NE.  N)  CO  TO  S 

JSU  •  JSUMINI 

DO  12  J  ■  l.JSU 

TARS  «  T  I  ISUM. J )  ♦  460.0 

Tf  *  Til  SUM, J I 

CALL  CONVITT,  XXHC,  EMISS,  XK I 

XXHR  •  EMISS  *  SIGMA  *  (TABS  ••  1  ♦  TABS  *•  2  *  ITSAB)  ♦ 

1  TABS  •  ITSAB)  **  2  ♦  ITSAB)  •*  31 

XXH  •  XXHC  »  XXHR 
IF (NKONT  .NE.  01  CO  TO  6 
PRINT  201,  XK.  XXHR,  XXH 

201  FORMAT  1 5H  XK  «,  E15.T,  5X,  TH  XXHR  •,  E15.T,  5X,  6H  XXH  »,  E15.7I 
6  02  1 1  SUM,  J I  <  XXH  /  XK  *  (TI1SUM.JI  -  TS) 

I F (NKONT  .NE.  01  GO  TO  12 
PRINT  10.  ISUM,  J,  021 1  SUM, J ) 

10  FORMAT  (5X.9H.  ISUM  •  ,  I3,3X,3H  J«.I5,  IOX,  11H  02IM,J)  »  ,  E12.4) 
12  CONTINUE 
CO  TO  90 
8  JSU  •  JSUMIKR I 
JSUI  ■  JSUMIKR, II 
I  FI JSU  .EQ.  JSUI)  CO  TO  90 
I F I JSU  .LI.  JSUI)  GO  TO  80 
JSUI1  «  JSUI  *  1 
DO  70  J  «  JSUI1,  JSU 
TABS  •  TIISU.J)  ,  460.0 
TT  >  TIISU,  Jl 

CALL  CONVITT,  XXHC.  EMISS,  XK I 

XXHR  •  EMISS  •  SIGMA  •  I  TABS  •*  3  »  TABS  ••  2  •  ITSAB)  » 

1  TABS  *  ITSAB)  ••  2  ,  ITSAB)  ••  3) 

XXH  «  XXHC  »  XXHR 

IF  I NKONT  .NE.  0)  CO  TO  9 

PRINT  201,  XK,  XXHR,  XXH 

9  021  ISU,  J)  •  XXH  /  XK  «  I  TIISU,  J)  -  TS) 

IF  (NKONT  .NE.  0)  CO  TO  70 
PRINT  203,  ISU.  J,  021 ISU, J) 

203  FORMAT (6H  ISU  •,  110,  SX,  4H  J  •,  1 10,  5X,  I2H  02IISU.JI  •.  E13.7I 
70  CONTINUE 
GO  TO  90 

80  JSUM1  •  JSU  ♦  1 
ISU1  •  ISU  ,  1 
00  85  J  •  J  SUM l ,  JSUI 
TABS  -  TIISU1.J)  *  460.0 
TT  •  TIISU1.J) 
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CALL  CONVITT,  XXHC,  EMISS 
XX HR  •  EMISS  *  SIGMA 
1  TABS  *  ITSAB) 

XXH  •  XXHC  ♦  XXHft 
IFINKONT  .ME.  0)  GO  TO  BS 
PRINT  201.  XK 
OZI ISU1.JI  *  XXH  /  XK 
69  CONTINUE 
90  CONTINUE 
RETURN 
ENO 


.  XXI 

•  (TABS  **  3  ♦  TABS  •*  2  •  ITSAB) 
••  2  ♦  ITSAB)  ••  J) 


*  ITS  -  Til SU1 • J) ) 


♦ 


SUBROUTINE  ORSUBI T.OR. TSAB. JM.LLS I 

COMMON  /BLKI/JSUMI8I,  N.  NKONTt  NP.  ISUM.  LI  I .  M,  SIGMA,  TS 

DIMENSION  T I  AO. 30 1 .  OR!  AO, JOI .LLSI  8 1  •  XI20I,  YI20I 

DU  20  1*1.  ISUM 

TABS  *  T I  1 ,  1 )  ♦  460.0 

TT  «  TII.ll 

CALL  CONVITT, XXHC, EMISS, XK) 

XX HR  •  EMISS  *  SIGMA  •  I  TABS  •*  J  ♦  TABS  ••  2  •  ITSAB)  ♦ 

1  TABS  *  (TSAB)  ••  2  ♦  ITSAB)  *•  31 

XXH  •  XXHC  *  XXHR 

20  OR  I  I ,  I )  «  XXH  /  XK  •  (TII.ll  -  TS) 

LI  •  LL1 
IS  *  0 

00  200  IK  * l , N 
JSU  •  JSUMI IK  I 
II  «  IS  «1 
IS  *  IS  »  LLSI IK ) 

DO  110  I  •  II.  IS 
TABS  •  Til, JSU)  ♦  460.0 
TT  *  HI. JSU) 

CALL  CONVITT,  XXHC,  EMISS,  XK  I 

XXHR  «  EMISS  •  SIGMA  •  (TABS  ••  1  •  TABS  ••  2  •  ITSAB)  » 

1  TABS  *  (TSAB)  ••  2  ♦  ITSAB)  ••  3) 

XXH  •  XXHC  ♦  XXHR 

130  OR  I  I  ,  JSU  I*  XXH  /  XK  •  (Til,  JSU)-  TS) 

200  CONTINUE 
RETURN 
END 


SUBROUTINE  AXIOERII,D2,OTZ.DTI2,Q2,ILSI 

DIMENSION  T 140, 301 •  OT2(40,301,  OT22I40,30),  02(51,  0ZI40.30I, 
1  LLSIBI.  X I  20 1 ,  VI 201 

COMMON  /BLK l/JSUNI 8 1 ,  N,  NKONT,  NP,  ISUM,  LL1 ,  M,  SIGMA.  TS 
IS  *  0 

00  320  IK  >1,  N 
JSU  •  JSUMI IK) 

C  IMIK  .EQ.  N)  GO  TO  150 
C  I F I JSU  .LE.  JSUMI IK4 1 | )  GO  TO  150 
II  •  IS  «1 
IS  *  IS  ♦  LLSI  IX I 
150  00  320  JH.JSU 

IF!  1 1  .NE.  1)  GO  TO  185 
02 1 1 ,  J  1*0.0 
01211, J)  — 0211, J) 

OT 221 1, J )  •  1./IS./DZI 1 1  ••  2  *1-85.  »  Tll.JI  »  108.0  •  TI2.J) 
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1  27.0  •  Tll.J)  ♦  4.0  *  TI4.JI1  ♦  11.  7  1.0  •  07(1. J)  7  *211) 

IFIMOOINKONT.NPI  .NE.  01  GO  TO  192 
PAINT  20.  J,  DT2U.JI,  DT22I1.J) 

20  FORMAT ( 7H  UTZ(1,,I1,  )H)  E19.6,  10X.  7M  6T|2  >.119. 6) 

192  12  *  II  1 

0T7II2.  1./6. 707(1X1  *  I  -2.  *  TII2-I.JI  -  3.  *  TII2.J)  »  6.  *T( 

112, J)  i ' 1 2*2, J 1 1 

0T22II2,JI  •  1.  702 (IK  I  ••  2  •  ITII2M.JI  -  2.  •  T 1 1 2 , J I  ♦  TII2-1, 
1  J  I  > 

IFIMOOINKONT.NPI  .NE.  0)  GO  TO  31 
PRINT  22,  12.  J.  0T2 1 1 2, J I .  DT72II2.J1 
22  FORMAT  1 9M  OT2I.I3.2H  ,  13,  3HI  •,  E19.6,  10X,  TH  OT22  *,  E19.A) 

31  I NT  *  II  *2 
IS,’  >  IS  -  2 
00  193  I  •  INT,  I S2 

DTZll.JI  »  1.  7  12.  7  D2IIKI*  (Til  -  2.  Jl  -  6.  *  Til  -  1,  Jl  ♦  8. 
1*  Til  ♦  1,  J)  -  Til  ♦  2,  J 1 1 

C  193  Of  22 1  I , J I  «  1.  7  12.  /  02  I  IK  )•»  2  •  I-TII-2.JI  *  16.  *  T( 1-1. Jl  - 

DT  22 1 I , J )  *  I.  7  12.  7  02(1X1**  2  •  I-TII-2.J)  ♦  16.  *  T( 1-1, Jl  - 

1  30.  •  TU.JI  ♦  16.  »  T(  1*1,  J I  -  T(  I  *2,  Jl  I 

IFIMOOINKONT.NPI  .NE.  01  GO  TO  193 

PPI  <  V.  I.  J.  0T2II.  Jl,  0TZ2II.  Jl 
I'J  CONTINUE 

111  =  IS2  ♦  1 

nrziui.Ji*  l.  /  6.  /  oziixi*  i  -6.  •  riin-i.Ji*  3.  •  tiiii.ji*  2 
i.*riiu*i,ji  *  r 1 1 1 1 - 2 . j 1 1 

01  72  I  I  I l . J  I  *  1.  7  07111  **  2  •  (T(I!1»1,J)  -  2.  •  T 1 1 1 1 , Jl *T 1 1 ll-l 
1,  Jl  I 

IFIMOOINKONT.NPI  .NE.  01  GO  TO  32 
PRINT  22.  II1.J.  DT7  I  I  1 1, J 1.01721  I  1 1, Jl 
12  IFIIK  .EC.  N)  GO  TO  1&9 

I F  I J  .GT.  JSUM I  I X  *  1 1  I  GO  TO  169 

DT  7  I  I S , J  I  ■  1./07 I |K*1)/  I  1 • *07  I  IK* 1 1 /  D2 I  I  X  1  I  *  I  T 1 1 S» 1 ,  Jl  -  Til 
1S.JI  •  I  1.-07  1 1 X ♦ 1 1**  2  7  07(1X1**  23-02  I  IK* 1 1 ••  2  7  OZIIKI**  2 

2  •  I  I  IS-1,  Jl  1 

OT72I1S.JI  •  2.  / 1 07 1  IK  I *07 1 IX* 1 1 *02 1 1 X* 1 1 **2 1 *l#2 1 1 X* 1 1 762 1 IX) • 

1  T  US-1,  JI-ID2  I  IX*l  17  02  I  IK  I  •  1.1  *  TIIS.JI  ♦  T 1 1  S*  1  .  J I  I 
IFIMOOINKONT.NPI  .NE.  0)  GO  TO  120 

PRINT  22.  IS.  J.  DT7IIS.J),  OT22IIS.J) 

GO  TO  120 

160  IFIIK  .EO.  1)  GO  TO  120 

I F I J  .IE.  JSUMIIX-lll  GO  TO  169 
169  0T7IIS.JI  *-07 I  I S, J I 

0T22IIS.J)  «  1.  /  18.  7  07(1X1**  2  •  (-19.  •  TIIS.JI  »  101.  •  TITS 

I- 1. Jl  -  27.  *  T I  IS- 2. J I  *  A.  •  TIIS-3.JTI  -  11.  •  07IIS.J)  7  3.  7 

2  07  UK  I 

IFIMOOINKONT.NPI  .NE.  01  GO  TO  320 
PRINT  22.  IS.  J.  0T7IIS.JI.  OT72I1S.JI 
GO  TO  320 

189  IFIJSU  .GT.  JSUMIIK-ll)  GO  TO  200 

190  DT7UI.J)  •  1.  /  6.  7  02(1X1*  (-2.  •  Tlll-I.J)  -  1.  •  TIll.JI*  6.* 

1  T  I  II  *  1  •  J I  -  T  II  I  *2, J I  I 

OT77III.J)  *  1.70212)  ••  2  *  IT(li*l,J)  -  2.  •  TIll.JI  *TIII-I,J)> 
IFIMOOINKONT.NPI  .NE.  01  GO  TO  33 
PRINT  22,  II,  J.  0T2I1I.JI,  0T22I  I  I . Jl 
M  GO  10  192 

200  IFIJ  .IE.  JSUMIIK-ll)  CO  TO  190 
0T7I I  I, J)  •  021  1 1 , J I 

0T22 1  1 1  *  J  I  •  1.7  18. 702  I  IK  I  •»  2  *1-89.  •  TIll.JI  »  108.0  *  TIIU1.J 

II- 27.0  •  Till *2, J I  *  4.0  •  Till *3, J 1 1  -  11.  •  3.0  •  02(11. J)  7  021 

2  IK) 
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1 F I  MODI NKONT  ,  NP I  .HE.  01  CO  TO  34 
PRINT  22,  II,  J,  OT  Z  1 1 1 ,  J  I ,  0T22(II,JI 
34  I  Ft  IK  .NE.  N)  GO  TO  140 
CO  TO  142 
320  CONTINUE 
'  RETURN 
END 


SUBROUTINE  R ADDER IT,DR,0TR,DTR2,0R,JJS,LLSt 

DIMENSION  T 1 40, 30),  DRI4I,  DTR(40,30I,  DTR2140.J0I,  QRI40.30I, 

I  JJSI8),  LLSI8I,  XI 20) ,  VI20I 

COMMON  /BLKl/JSUMIBI,  N.  NKONT,  NP,  I  SUM,  LL1,  M,  SIGMA,  TS 
C 

XH  *  200.0 
TG  •  2000.0 
C 

1 1  r  *  i 
I  LA  •  LLSIi) 

DU  360  IK  •  1,  N 
JSX  «  JSUMIIK) 

00  104  JR  •  l,N 
JSX  *  JSX  -  JJSIJR) 

IF  I JSX  .EO.  01  GO  TO  306 
304  CONTINUE 
JR  «  M 

306  CONTINUE 

DO  340  I  •  IFF, I L  A 
TT  *  TII.l) 

CALL  CONVITT.  XXHC,  EMISS,  XKI 
OR  I  I  •  1  l«XH  *  I  TG  -  Til. 1)1  /  XK 
OTRII.lt  — QRII.ll 
C  ORI  1,1)  «  7  7 

0TR2II.I)  «  1.  /  18.  /  OR  I  1 )  •*  2  *  (-84.  •  TII.l)  4  108.  •  Til, 2) 
1  -  27.  *  Til, 31  ♦  4.  *  TII.4))  ♦  II.  /  3.  /  DRI1)  *  QRII.ll 
I F I  MOO  I  NKONT,  NP  I  .NE.  01  GO  TO  40 
PRINT  20,  I,  OTRI  1,1),'  DTR2I  1,1) 

20  FORMAT  1 4H  OTRI, 13.  4H.  1 )  E14.6.  10X.  7H  0TR2  *,  E14.6) 

40  OTRI  1,21  *  1.  /  6.  /  OAU)  •  1-2.  •  TII.l)  -  3.  *  Til, 21  ♦  6.  • 

1  Til, 3)  -  TII.4)) 

0TR2II.2)  ■  1.  /  OR  I  1 )  ••  2  •  ITU, 3)  -  2.  •  Til, 2)  ♦  Til, II) 
IFIMOOINKONT.NP)  .NE.  0)  GO  TO  41 
PRINT  22,  1,  OTRI 1 . 2  1 ,  0TR2I1.2) 

22  FORMAT  1 4H  OTRI, II,  4H,2)  •,  E15.6,  10X,  7H  DTR2  «,  E14.6) 

41  JJ1  ■  JJSI1I  -  2 
00  310  J  •  3,  JJ1 

OTR(I.J)  •  1.  /  12.  /  DRI 1 1  *  I T I  I . J— 2 >  -  8.  •  TII.J-l)  4  8.* 

1  TII.J41)  -  T I  I , J*  21  I 

0TR2II.JI  «  1.  /  12.  /  DRI1I  **  2  •  l-TII.J-2)  ♦  16.  •  TII.J-ll  - 
1  30.  •  TII.J)  »  16.  »  TII,J41)  -  TII.J42)) 

IFIMOOINKONT.NP)  .NE.  0)  GO  TO  310 
PRINT  24,  I,  J,  DTRII.J).  DTR2II.J) 

310  CONTINUE 

24  FORMAT  1 5H  OTRI.13.2H  ,13, 5H)  «,  E14.6,  10X,  7H  DTR2  «,  E14.6) 

0TRII.JJ141)  •  1.  /  6.  /  DRI 1 1  *  1-6.  *  TII.JJl)  4  3.  *  T 1 1 , J J 1 4 1 ) 
1  42.  •  TII.JJ142)  4  T 1 1 , J J 1— 1)1 

0TR2 1  I , J J 14 1 1  ■  1.  /  ORI 1 1  ••  2  •  ITII.JJ142)  -2.  *  Tl I , JJ141I4TI I 
1.  JJ1 ) ) 

IFIMOOINKONT.NP)  .NE.  01  GO  TO  43 
JZ  •  JJl  41 
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PRINT  24,  I,  JZ,  DTR 1 1 , JZ ) ,  0TR2I1.JZ) 

43  IFIJJSIl)  .60.  JSUMIIKI)  GO  TO  325 
JJ  >  JJSIl) 

DTRII.JJ)  •  1.  /  OR  1 2 1  /  II.  ♦  OR  1 2)  /  DRIll)  •  (T(1,JJ*1>  -  Til, 

1  JJ)  *11.  -  OR  1 2 1  ••  2  /  ORI 1 )  ••  2)  -  0RI2)  ••  2  /  DRI1I  *•  2  * 

2  T 1 1 , J J- 1 1  ) 

0TR2II.JJ)  »  2.  /  I  ORI 1 )  •  ORI 2)  ♦  0RI2I  ••  2)  •  IDRI2I  /  ORIll  * 

1  Til, JJ- 1 1  -  I ORI 2 )  /  ORIll  »  II  *  Til , JJ)  »  TII.JJ  *111 
IF (MODI NKONT.NP I  .NE.  01  GO  TO  44 
PRINT  24,  I,  JJ,  OTRII.JJI,  0TR2II.JJI 

44  JSX  *  0 

00  315  JK  •  2,  JR 
JSX  *  JSX  ♦  J JS I JK*  1 1 
JP  «  JSX  ♦  l 

DTRII.JPI  *  1.  /  6.  /  ORI  JK)*  1-2.  •  T(  I  ,  JP- 11-3.  •  TII.JP)  ♦  6.  * 

I  T 1 1 , JP* 1 )  -  Tl  I,JP*2I) 

0TR2II.JP)  >  1.  /  OR  I JX  )••  2  •  ( Tl I , JP* 1 1  -  2.  •  TII.JP)  »TII,JP-t 

II  ) 

IFIHOOI NKONT.NP)  .NE.  0)  GO  TO  45 
PRINT  24,  I,  JP,  OTRII.JPI,  0TR2I I , JP) 

45  Jl  *  JP  ♦  l 

JL  «  Jl  ♦  J JS I JK  I  -  4 
00  312  J  •  JI.JL 

OTRII.J)  «  1.  /  12.  /  ORI JK 1  *  ITII.J-2)  -  S.  *  TII.J-l)  »  8.  * 

1  Tl  1,  J*1 1  -  Tl  I,  J  *  2  >  ) 

DTR2II.J)  «  1.  /  12.  /  OR  I JK )**  2  *  l-TII.J-2)  »  16.  •  T 1 1  • J- 1 1  - 
1  30.  •  T I  I . J  )  »  16.  *  T 1  I , J*1 1  -  Tl I , J*2) ) 

I F I NOOI NKONT , NP )  .NE.  0>  GO  TO  312 
PR1NI  24,  I,  J,  OTRII.J),  0TR2II.JI 
312  CONTINUE 

JPP  -  JL  ♦  1 

OTRII.JPPI  -  1./  6.  /  DRIJK1*  (-6.  *  Til, JPP- 1 1  ♦  3.  •  Til, JPP)  ♦ 

1  2.  *  T I  I , JPP* 1 )  ♦  Til, JPP-2 1 ) 

0TR2II.JPP)  »  1.  /  DR  I JK I**  2  «  I  Tl  I , JPP*1 )  -  2.  *  TII.JPP)  ♦  Til, 
1  JPP- 313 

IFIHOOINKONT.NPI  .NE.  0)  GO  TO  47 
PRINT  24,  I,  JPP.OTRI  I  ,JPP),0TR2I  I, JPP) 

47  IFIJK  .EO.  JR)  GO  TO  315 
JPL  «  JPP  ♦  1 
JZ  «  JK  ♦  l 

OTRII.JPLI-  1.  /  OR  I JZ  )  /  11.  ♦  ORI  JZ )/  ORIJKX*  IT(I.JPL*1I-  Til, 

1  JPL  )  *  I  1 .  -  OR  I  JZ  )  *•  2  /  ORI  JK  )•*  21  -  ORIJZ)**  2  /  ORIJK)**  2  * 

2  TII.JPL-D) 

0TR2II.JPLI-  2.  /  IDRIJK)*  ORIJZI*  ORIJZ)**  21  •  IDRIJZ)/  ORIJK)* 

1  TII.JPL-ll-  I  OR  I  JZ  I  /  ORI  JK  I  ♦  II  •  Tlt.JPU*  TII.JPL*!)) 

I F IHODINKONT , NP )  -NE.  0)  GO  TO  315 
PRINT  24,  I,  JPL.OTRI  l,JPLI,0TR2l I.JPL)  ■ 

315  CONTINUE 

325  JLT  *  JSUMtIK)  -  2 

DTR 1 1 , JL  T  *2  I  *-QRI l,JLT*2) 

0TR2I  I,  JLT*2)  •  1.  /  18.  /OR  I  JR )  *•  2  *  1-85.  *  T1I,JLT*2I  *  108.  • 

1  Til,  JL  T*  l )  -  27.  *  Til,  JLT)  ♦  4.  •  TII.JLT-D)  -  11.  2  3./  ORI 
2 JR ) *  OR  1 1 ,  JLT*2I 
IFIHOOINKONT.NPI  .NE.  0)  GO  TO  350 
JW  *  JLT  ♦  2 

PRINT  24,  I,  JM,  OTRII.JUI,  DTR2I 1 , JW) 

350  CONTINUE 

1FIIK  .EQ.  N)  GO  TO  360 
IUN  «  IK  *  1 
I  FT  •  1FT  ♦  LLSIIK) 

ILA  •  IL*  *  LLSI IUN) 


28 


360  CONTINUE 
RETURN 
END 


SUBROUTINE  DKOT I  TT»  XK»DKT, KRHOi XCR  »  X.  VI 
DIMENSION  X 1 20) •  VI20) 

CALI  XKKSI TT,  XK , DKT ) 

XRHO  *  490.0 
I  ■  1 

CALL  LINEARITT,X,Y«  XCP,  I ) 

RETURN 

END 


SUBROUTINE  XKKSI TT. XX. OKT) 
IF1TT  .GT.  1472.01  GO  TO  14 
XK«28.30-.00870*TT 
DKT  •  -  0.0087 
GO  TO  20 

14  XK«10.39».00347*TT 
DKT  •  0.00347 

20  CONTINUE 
RETURN 
ENO 


SUBROUTINE  L INEARI A, X, V, VV,  1 1 
DIMENSION  X ( 20 1 >  Y( 20) * 

1  IFIYIIMI  .IT.  V(l)l  BO  TO  too 

C  USE  FOLLOWING  IE  AS  V  INCREASES  X  INCREASES 

10  1FIA-XIIII3.2.2 

C  USE  FOLLOWING  IF  AS  V  INCREASES  X  DECREASES 

100  IFIA-XI I ) 12,2,3 

2  1-1*1 
GO  TO  1 

3  I-I-l 

VV-YI  I  l*(A-XC 1*1)  l/IXIII-XI 1*1) »*Y«  l«ll*(A-X(I ) )/<X(l«l)-X(l>) 

RETURN 

ENO 


3 

4 

1714E-G8  70.0  20.0 

99550 

200 

7 

5 

5 

5 

5 

5  8 

.052083 

.0703125 

.078125 

.09896 

.26333 

1.0417 

1.81907 

3.50 

3 

2 

l  1 

0.0 

.108 

200.0 

.112  400.0 

.125 

600.0 

.132 

800.0 

.150 

1000.0 

.160  1200.0 

.185 

1600.0 

.180 

2000.0 

.180 

2200.0 

.150 
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I 


I 


M  a 

JJS 


SIGMA  • 

4 

4 


n.iriAOOf-09 

» 

4 


fS  •  O.TOOOGE  <9? 


7f|Mf  •  o.^ogooe  #>7  • ' 


O.WOI  lOOf-ol  O.JOlMiOf-OI  0./.W(00(-01  l.Nmt.,1 
0.J6JJJ00C  00  0.|P«l699(  01  0.1(19070(01  0.7(00000(  01 

f 


KM 
JSUM(l) 
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TIME 


0.2  T06E-02 


A*lAl  LOCATION 
0.76ME  0?  n 

0.7000E  02  0 

0.7000E  0? 
axial  location  • 

0.76ME  02  0 

O.TOOOE  0?  0, 

0.7C00F  02 
AXIAL  LOCATION  - 
0. 76*  IE  02  0, 

O.TOOOE  02  0. 

0.7000E  02 
AXIAL  LOCATION  - 
O.TMir  02  0, 

O.llMMlf  l»2  |l, 

o,/mini  *i2 

AKIAl  LPCAIIIIN  . 

0.76AIE  02  r. 

0 «  7000E  C2  0. 

O.TOOOE  02 
AXIAL  LUflATION  - 
0.7#.*  II  02  it. 

(I.  loom  02  ,  II, 

A4IAL  UK  A|  |i,n  - 
O.K.AH  02  0. 

0. 7000E  02  0. 

AXIAL  LOCATION  - 
0. 76*  IE  02  o, 

0.7000F  02  n. 

AXIAL  LUCATIOn  - 
0.76AIE  02  r. 

0.7000E  02  0. 

AXIAL  L0CA1  IriN  - 
0.76ME  02  0. 

0.7000*  02  0. 

axial  location  - 

0.76*  IE  0*  C, 

AXIAL  LOCATION  - 
0.7ft*lE  C2 
AXIAL  LIICATIIIN  - 
0. 76*  |  f  02  0. 

AXIAL  LOCATION  - 
0 ( 76* 1 E  02  i 

axial  location 


7000E  62 
700CE  02 


.  7000C  02 
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FLOW  CHART  OF 
MAIN  PROGRAM 


LIST  OF  SYMBOLS  USED  IN  TEXT 
Symbols 

C  -  specific  heat  (BTU/lb  *F) 
h  -  heat  transfer  coefficient  (BTU/hr  ft2  *F) 

1  -  number  of  radial  nodes 

j  -  number  of  axial  nodes 

k  -  thermal  conductivity  (BTU/hr  ft  °F) 

m  -  axial  node 

n  -  radial  node 

q"  -  heat  flux  (BTU/hr  ft2) 

R  -  radial  boundary  (ft) 
r  -  radial  coordinate  (ft) 

T  -  temperature  (*F) 
t  -  time  (hr) 

Z  -  axial  coordinate  (ft) 
aR  -  radial  Increment  (ft) 
a Z  -  axial  Increment  (ft) 
e  -  emlsslvlty 
p  -  density  (lb/ft3) 

o  -  radiation  coefficient  (0.1714  BTU/hr  ft2  °R‘‘) 
Subscripts 
9  -  9«s 

1  -  Initial  value 
m  -  node  m 

n  -  node  n 

0  -  surroundings 

*  •  boundary  1 ,  sequent  1 

2  -  boundary  2,  segment  2 

3  -  boundary  3,  segment  3 
w  -  boundary  4 
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